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INTRODUCTION 


A  perforated  muzzle  brake  consists  simply  of  a  set  of  vent  holes  drilled 
through  the  wall  of  a  cannon  near  the  muzzle  (see  Figure  1).  Compared  with  con¬ 
ventional  baffle  brakes,  they  are  lighter  and  simpler  to  manufacture  and,  as 
shown  in  a  series  of  reports  by  Dillon  and  Nagamatsu  (refs  1-5),  they  can  be 
designed  to  provide  significant  levels  of  recoil  reduction.  Also,  because  the 
vented  area  can  be  located  symmetrically  around  the  tube,  a  more  favorable  flow 
environment  is  provided  for  finned  projectiles.  This  is  an  important  con¬ 
sideration  for  weapon  accuracy  and  the  structural  integrity  of  the  projectile. 

In  a  field  study  of  105-mm  brake  designs,  it  was  found  that  asymmetrical  venting 
can  lead  to  bending  and  even  breakage  of  the  fins. 


Figure  1,  Schematic  drawing  of  a  perforated  muzzle  brake. 

Nagamatsu,  Duffy,  Choi,  and  Carofano  (ref  6)  made  a  numerical  calculation 
of  the  steady  three-dimensional  flow  through  a  vent  hole  in  the  wall  of  a  shock 
tube.  The  predicted  pressure  distribution  on  the  vent  wall  compared  favorably 
with  the  experimental  measurements  of  Nagamatsu,  Duffy,  and  Choi  (ref  7).  It 
was  also  shown  that  these  results  could  be  combined  with  a  one-dimensional  model 
of  the  transient  flow  in  a  cannon  to  predict  the  impulse  reduction  produced  by  a 
perforated  muzzle  brake.  The  predictions  agreed  well  with  the  experimental 


measurements  of  Dillon  (ref  1)  for  a  20-mm  cannon. 


These  calculations  were  made  for  a  perfect  gas,  but  at  the  pressure  levels 
found  in  larger  caliber  weapons,  the  necessity  of  including  the  gas  covolume 
correction  in  internal  ballistics  solutions  is  well  known  (ref  8).  The  impor¬ 
tance  of  this  factor  on  perforated  muzzle  brake  performance  is  considered  here. 

It  is  also  of  interest  to  know  if  and  how  experimental  results  obtained 
with  a  small  caliber  weapon  can  be  scaled  up  in  an  effort  to  avoid  the  high 
costs  associated  with  testing  larger  weapons.  In  studies  of  conventional  muzzle 
brakes  (refs  9-11),  the  gasdynamic  efficiency,  /5,  has  been  found  useful  for  this 
purpose.  Its  application  to  perforated  muzzle  brakes  is  discussed  here. 

THE  THREE-DIMENSIONAL  MODEL 

When  the  propellant  gas  expands  through  the  brake,  an  asymmetric  pressure 
distribution  develops  in  each  hole  with  the  highest  pressures  acting  on  the 
downstream  surface.  The  vector  and  pressure  contour  plots  of  Figures  2  and  3 
show  typical  flow  patterns  in  the  symmetry  plane  of  one  hole  and  the  portion  of 
the  tube  associated  with  it.  The  flow  variables  in  the  tube  are  uniform  across 
the  entrance  plane.  The  solid  lines  in  the  vector  plot  indicate  where  the  local 
Mach  number  is  unity. 

In  Figure  2,  the  flow  enters  at  Mach  one  and  accelerates  to  supersonic 
velocities  as  a  portion  of  the  gas  expands  and  turns  into  the  hole.  The  shock 
at  the  downstream  lip  of  the  hole  turns  the  expanded  flow  parallel  to  the  solid 
surfaces  and  reduces  the  velocities  to  subsonic  levels.  The  pressure  on  the  lip 
is  nearly  twice  the  static  pressure  of  the  incoming  stream.  The  flow  acceler¬ 
ates  away  from  this  region  and  leaves  the  tube  and  hole  at  supersonic  veloci¬ 
ties.  There  is  a  large  subsonic  region  on  the  upstream  portion  of  the  hole 
where  the  pressure  is  nearly  uniform.  A  more  complete  description  of  the  three- 


Fortunately,  the  flow  contains  many  features  which  permit  a  vigorous  simplifica¬ 
tion  of  the  problem. 


Figure  3.  Velocity  vector  and  pressure  contour  plots  for  a  flow  with  a  Mach 
number  of  two  at  the  entrance  plane. 

i  First,  because  of  the  large  volume  of  the  gun  tube,  the  blowdown  process 

I 

takes  on  the  order  of  tens  of  milliseconds  while  the  three-dimensional  calcula¬ 
tions  indicate  that  the  flow  in  a  hole  is  established  in  a  fraction  of  a  milli- 

j  second.  Therefore,  the  latter  can  be  treated  as  quasi-steady  and  only  the  flow 

1 

'j  within  the  tube  must  be  considered  as  a  time-dependent  process. 

Secondly,  in  the  applications  of  interest,  the  flow  is  either  sonic  or 
supersonic  as  it  enters  the  brake  and,  due  to  the  venting,  expands  to  higher 
Mach  numbers  as  it  travels  downstream.  Also,  because  of  the  high  tube 


4 


pressures,  the  gas  exits  each  hole  near  sonic  or  supersonic  velocity  over  most 
of  the  exit  plane  area  (see  Figures  2  and  3).  Experience  has  shown  that  the 
flow  is  rather  insensitive  to  the  outflow  boundary  condition  over  the  remaining 
subsonic  portion.  Thus,  the  flow  at  a  particular  hole  location  is  not 
influenced  by  events  occurring  farther  downstream  or  outside  of  the  tube.  It 
depends  solely  on  the  conditions  in  the  tube  upstream  of  the  hole.  This  obser¬ 
vation  permits  a  dramatic  reduction  in  the  amount  of  three-dimensional  com¬ 
putation  required  to  predict  brake  performance. 

The  Euler  equations  may  be  written  in  conservative  form  as 

aq  3F  dG  dH 

—  +  —  +  —  +  —  =  0  (1) 
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p  is  the  density;  m  =  pu,  n  =  pv,  and  J  =  pw  are  the  momentum  components  in  the 
X,  y  and  z  directions,  respectively;  u,  v,  and  w  are  the  corresponding  velocity 
components.  P  is  the  pressure  and  E  is  the  total  energy  per  unit  volume  defined 
as 


E  =  pe  (m*+n*  +  l*)/2p  (2) 

where  e  is  the  specific  internal  energy. 

To  include  the  influence  of  the  gas  covolume,  n»  the  Abel  equation  of  state 
is  used.  For  this  model,  the  pressure  and  soundspeed,  c,  are  related  to  the 
state  variables  p  and  e  by  the  expressions 
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p  =  (Y-l)pe/{l-np) 
c*  =  y(y-l)e/(i-r)p)* 


(3) 


(4) 

where  y  is  the  specific  heat  ratio. 

Consider  the  result  of  nondimensional izing  the  Euler  equations  in  the 
following  way; 

p'  =  p/P2  ,  m'  =  m/y'P2P2  .  n'  =  n/yP2P2  »  =  l/yP2P2 

P'  =  P/P2  ,  E'  =  E/P2  ,  e'  =  ep2/P2  -  C  =  C/P2/P2 
x'  =  x/D  ,  y'  =  y/D  ,  z'  =  z/D  ,  f  =  t/p^7^/D 
where  P2  and  P2  are  the  density  and  pressure  of  the  uniform  flow  at  the  upstream 
plane  of  the  tube  and  D  is  the  vent  diameter.  The  form  of  the  Euler  equations 
remains  unchanged  while  the  inflow  boundary  conditions  become 


P'  =  1  (5) 

m'  =  M2yy/(l-nP2)  (6) 

n’  =  1'  =  0  (7) 

E'  =  (l-nP2)/(l-y)  +  yM2*/2(i-nP2)  (8) 

The  pressure  and  soundspeed  relations  become 

P'  =  (y-l)p’e'/(l-nP2P' )  (9) 

C*  =  y(y-l)e'/(l-nP2P’ )*  (lO) 


Since  the  flow  depends  only  on  the  inflow  boundary  conditions,  which  are 
seen  to  be  completely  described  by  the  upstream  Mach  number,  M2,  the  specific 
heat  ratio,  y,  and  the  dimensionless  product,  r7p2,  one  solution  with  these 
parameters  specified  is  valid  for  all  upstream  pressures  and  densities. 

Although  a  wide  range  of  physical  states  is  encountered  during  blowdown,  it  is 
found  that  only  a  few  three-dimensional  solutions  are  required  to  describe  them. 

Marten's  Total  Variation  Diminishing  scheme  (ref  12)  was  used  in  conjunc¬ 
tion  with  a  time-splitting  algorithm  to  solve  the  Euler  equations  after  they 


were  transformed  to  a  generalized  grid.  Separate  grids  were  used  for  the  vent 
and  the  tube.  A  cylindrical  grid  was  used  for  the  vent  with  a  Cartesian  gnd  at 
its  center  to  avoid  the  geometrical  singularity  there.  The  circular  tube  was 
replaced  with  a  rectangular  channel  to  simplify  interpolation  between  the  two 
grids.  As  can  be  seen  in  the  vector  plot  of  Figure  2,  19  grid  points  were  used 
across  the  vent  diameter  and  for  the  azimuthal  coordinate.  The  length  of  the 
rectangular  channel  (the  primary  flow  direction)  was  represented  with  28  grid 

points,  its  height  by  24  points,  and  its  depth  (into  the  paper)  by  13  points. 

Ct  course,  the  size  of  the  channel  could  affect  the  flow  through  the  vent. 
To  check  this  possibility,  the  example  shown  in  Figure  2  was  repeated  with  the 
number  of  grid  points  used  to  represent  the  height  and  depth  of  the  channel 
approximately  doubled.  The  effect  on  the  solution  in  the  vent  was  minimal, 
indicating  that  the  channel  was  essentially  an  infinite  reservoir.  The  inflow 
Mach  number  in  this  example  was  unity.  At  Mach  two,  the  effect  would  be 

expected  to  be  even  less  because  the  disturbance  produced  by  the  vent  occupies 

much  less  of  the  channel  (see  Figure  3). 


THE  ONE-DIMENSIONAL  MODEL 

The  one-dimensional  equations  of  motion  with  venting  at  the  tube  wall  are 

3p  3m  1  dm 

* -  (11) 

3t  3x  A  dx 


3m  3(m*/p+P) 

m  dm 

3t  3x 

pA  dx 

3E 

3(m(E-^P)/p) 

(E+P)  dm 

3t 

3x 

pA  dx 
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The  vent  term  (l/A)dm/dx  represents  the  mass  of  fluid  per  unit  time  per  unit 
volume  leaving  the  tube  at  x.  A  is  the  bore  area.  The  fluid  is  assumed  to 
leave  at  the  local  velocity  u  in  the  momentum  equation  and  with  the  local 
enthalpy  per  unit  mass  {E+P)/p  in  the  energy  equation  (see  Reference  13  for  the 
derivation  of  these  equations). 

Data  from  the  three-dimensional  solution  are  used  to  obtain  an  average 
value  for  the  dimensionless  mass  flux  leaving  the  vent  using  the  expression 

=  (IAh)/*  p'w’dA  (14) 

where  w'  is  the  dimensionless  velocity  in  the  z-direction,  i.e.,  parallel  to  the 
vent  axis.  The  integration  is  carried  out  over  the  vent  exit  area,  A^.  The 
flux  is  applied  instantaneously  at  each  axial  location  in  the  one-dimensional 
model.  It  is  a  function  of  the  three  parameters  that  appear  in  the  three- 
dimensional  solution. 

The  mass  flux  also  depends  on  the  aspect  ratio  of  the  hole,  which  is 
defined  as  the  ratio  of  its  height,  L,  to  its  diameter,  D.  In  Reference  6,  the 
brake  configurations  were  limited  to  sets  of  uniform  diameter  holes.  However, 
from  the  structural  point  of  view,  it  may  be  advantageous  to  use  progressively 
smaller  holes  away  from  the  muzzle,  since  the  web  of  material  between  each  suc¬ 
cessive  row  of  holes  must  support  the  full  load  generated  by  all  of  the  holes 
farther  downstream.  To  allow  for  this  possibility,  the  h^le  diameter  will  be 
taken  as  the  following  function  of  distance,  y,  along  the  brake: 

^  “  ^min  (Omax'^min)  0<a<l  ,  0<y<L\/  (15) 

where  Ly  is  the  length  of  the  vented  region,  Dmin  is  the  diameter  of  the  first 
row  of  holes  at  the  brake  entrance,  and  is  the  diameter  of  the  last  row  of 

holes  near  the  muzzle.  Let  Np  be  the  number  of  rows  of  holes  (rows  run  around 


the  circumference  of  the  tube)  and  Nc  be  the  number  of  columns  of  holes  (columns 
run  parallel  to  the  tube  axis).  For  a  uniform  spacing  of  the  holes  along  the 
brake,  the  vent  area  per  unit  length  is  NQNRnO*/4Lv  and  the  total  vent  area.  Ay, 
is  found  by  integration  to  be 

Ay  =  (7TNcNRDminV4)(l+2r/{a+l)  +  r*/(2a+l))  (16) 

where  r  is  the  ratio 

^  “  (^max~^min)/^min 

The  vent  term  can  be  written  in  dimensional  form  as 

e 

(l/A)dm/dx  X  -p'w'/Pp  (NqNr/Lv) (D/Db)*  (18) 

where  Db  is  the  tube  diameter.  P  and  p  are  the  local  values  of  pressure  and 
density  in  the  one-dimensional  solution;  they  appear  since  these  quantities  were 
used  to  nondimensional ize  the  product  p'w*  in  the  three-dimensional  solution. 

Using  a  control  volume  drawn  around  the  outer  boundaries  of  the  tube, 
including  the  muzzle  exit  plane,  the  following  equation  for  the  impulse,  I,  is 
obtained: 

I  X  WVpe  *  a/*  (Pe-Po+PeUeMdt  ♦  a/"[NcNr/Lv)/^'^p’w’u-  P(D/Db)  *dy]dt  (19) 
The  first  term  is  the  projectile  momentum  as  it  passes  the  brake  exit  plane  with 
velocity  Vpg,  The  second  term  is  the  thrust  produced  by  the  gas  discharged 
through  this  plane.  The  subscript  "e"  refers  to  time-dependent  properties 
obtained  from  the  one-dimensional  solution.  Pq  is  atmospheric  pressure. 

The  quantity  in  brackets  in  the  third  term  represents  the  axial  thrust  pro¬ 
duced  by  the  gas  venting  through  the  holes.  As  can  be  seen  from  the  vector 
plots  in  Figures  2  and  3,  the  gas  velocity  is  not,  in  general,  parallel  to  the 
hole  axis  everywhere  across  the  exit  plane,  so  this  integral  may  be  negative, 
zero,  or  even  positive.  The  average  value  of  the  momentum  flux  is  calculated 
from  the  three-dimensional  solution  using  the  expression 
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(20) 
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p'w'u'  s  (l/A^^)/^^p'w' u'dA 
where  u'  is  the  dimensionless  axial  velocity  component  of  the  vented  gas.  This 
quantity  is  a  function  of  the  local  Mach  number  and  geometry,  just  as  with  the 
mass  flux.  The  local  values  of  pressure  and  hole  diameter  appear  in  the  spatial 
integral  in  Eq.  (19)  because  they  are  position-dependent.  The  integral  is  eval¬ 
uated  at  each  time-step  in  the  one-dimensional  solution. 

The  calculation  is  started  with  the  projectile  base  just  upstream  of  the 
vented  region.  Initial  conditions  are  taken  from  an  internal  ballistics  code 
which  uses  the  Pidduck-Kent  limiting  solution  (ref  8)  to  compute  the  gradients 
of  pressure,  density,  and  velocity.  These  are  given  by 


Q  =  iiill  c  c  ^  (l:Qy!)l/(y-l)^  )-l 

2y  W  '  2y  W  1-Q 

X  =  x/xs 

z  -  n  §  (Z  -  (liQy!)i/(Y-i)d4) 

u  0  i-M  0 

p  *  (.l:9.)-Y/(y-i) 

V  li-QZ*' 

1  .  n  =  (U  -  n)  ?  (.l:?.)i/(r-i) 

p  '  'c  "  Zy  £2  W  'l-QZ*' 

u  =  ZVpv 


(21) 

(22) 

(23) 

(24) 

(25) 

(26) 


where  Py  is  the  projectile  base  pressure  at  the  starting  time,  ts.  and  Vpv  is 
the  projectile  velocity.  C  and  W  are  the  charge  and  projectile  masses,  respec¬ 
tively.  n  is  a  parameter  which  depends  on  the  ratio  C/W  and  y.  It  is  evaluated 
once  at  the  start  of  the  solution  by  solving  Eq.  (21)  by  iteration.  Xg  is  the 
projectile  position  at  time  tg;  it  includes  a  uniform  extension  of  the  tube  of 
length  U/A  where  U  is  the  chamber  volume  and  A  is  the  bore  area. 


i 


t 
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At  each  grid  point  in  the  one-dimensional  solution,  X  is  computed,  then  Z 
is  found  by  iteration  from  Eq.  (23).  The  values  of  pressure,  density,  and 
velocity  are  then  calculated  and  the  conservative  variables  in  the  Euler 
equations  are  formed  from  these. 

The  equation  of  motion  for  the  projectile  is  solved  together  with  the  con¬ 
servation  equations  until  the  projectile  exits  the  barrel.  The  effect  of  brake 
venting  on  the  final  projectile  velocity  is  thus  determined.  The  remainder  of 
the  calculation  deals  with  the  motion  of  the  gas  in  the  barrel  until  complete 
blowdown  is  achieved. 

THE  DIMENSIONLESS  FLUX  FUNCTIONS 

Four  factors  affect  the  values  of  the  dimensionless  flux  functions:  the 
specific  heat  ratio,  y;  the  covolume  parameter,  r)P2;  the  upstream  Mach  number, 
M2;  and  the  vent  geometry,  y  does  not  vary  much  for  different  propellants  and 
its  effect  on  the  three-dimensional  solution  has  been  found  through  experience 
to  be  small.  A  single  value  of  y  *  1.22  was  used  here  to  reduce  the  amount  of 
calculation. 

A  study  of  internal  ballistics  solutions  for  various  caliber  weapons  showed 
that  the  covolume  parameter  is  not  likely  to  exceed  0.2,  therefore,  this  study 
used  only  this  value  and  the  perfect  gas  limit,  nP2  =  0.0. 

In  an  earlier  study  (ref  6),  the  flux  functions  were  found  to  vary  nearly 
linearly  with  M2,  therefore,  the  calculations  were  limited  to  the  two  cases  M2  = 
1  and  M2  =  2. 

The  hole  geometry  is  characterized  by  the  ratio  of  its  height,  L,  to  its 
diameter,  0.  The  height  is  equal  to  the  tube  wall  thickness  which,  in  most 
designs,  would  be  uniform  throughout  the  brake.  Manufacturing  considerations 
suggest  a  range  of  L/D  from  1  to  2  so  these  two  values  were  used  here. 
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These  parameter  choices  required  a  total  of  eight  three-dimensional  solu¬ 


tions.  The  averaged  flux  functions  are  shown  in  Figure  4. 


X 

D 

-JO  - 

li  N 


L/D=l 


L/D=2 


1.00  1.50  2.00 

MflCH  NUMBER 


L/D=2 


L/D=l 


l.OO  l.SO-  2.00 
MRCH  NUMBER 


Figure  4.  Averaged  mass  and  momentum  flux  functions  computed 
from  the  three-dimensional  solutions. 

The  dashed  lines  represent  the  effect  of  including  the  covolume  term  in  the 
equation  of  state.  Since  the  covolume  parameter  was  chosen  to  represent  an 
upper  limit  of  conditions  expected  in  a  brake,  it  would  not  appear  to  be  very 
important.  This  was  confirmed  with  several  test  cases  using  the  one-dimensional 
code.  It  can  be  concluded  that  in  any  future  three-dimensional  studies,  the 
effect  of  covolume  on  the  flux  functions  and  brake  performance  is  too  small  to 
warrant  doubling  the  amount  of  computation.  The  covolume  term  should  be 
included  in  the  equation  of  state  used  in  the  one-dimensional  code,  however,  to 
be  consistent  with  the  internal  ballistics  code  used  to  generate  the  starting 
solution. 

Note  that  the  momentum  flux  is  negative  for  both  vent  heights,  especially 
the  shorter  one.  This  is  consistent  with  the  velocity  vector  plots  of  Figures  2 
and  3  —  the  flow  leaving  the  shorter  vent  is,  on  balance,  directed  more 
upstream. 


WiWiW.Wi' 


The  linear  functions  shown  in  Figure  4  were  incorporated  into  the  one¬ 
dimensional  code  as  simple  analytical  expressions.  Interpolation  was  used 


wherever  necessary.  The  covolume  correction  was  included  since  the  information 
had  already  been  obtained. 

BRAKE  PERFORMANCE  VERSUS  VENT  AREA  DISTRIBUTION 

Previous  work  (refs  1-6)  has  shown  that  when  a  brake  is  made  with  a  set  of 
uniform  diameter  holes,  the  hole  diameter,  D,  should  be  set  equal  to  the  tube 
wall  thickness,  L,  to  achieve  the  greatest  reduction  in  weapon  impulse  for  a 
given  vented  area.  This  is  a  result  of  the  asymmetrical  pressure  distribution 
in  the  hole.  With  reference  to  Figure  2,  as  the  gas  expands  to  supersonic 
velocities  along  the  downstream  wall,  the  gas  pressure  falls  below  the  level 
acting  on  the  adjacent  upstream  surface.  The  crossover  occurs  at  a  distance 
into  the  hole  of  about  one  diameter.  Beyond  this  point  the  net  braking  load 
decreases.  Thus,  the  most  efficient  brake  is  a  set  of  uniform  holes  with  L/D  = 
1.  However,  a  less  efficient  brake  might  be  preferred  or  even  required  based  on 
structural  considerations.  This  will  be  demonstrated  using  a  120-mm  cannon  as 
an  example.  The  values  in  Table  I  were  used  to  start  the  solution.  L[j  is  the 
barrel  length,  i.e.,  the  total  distance  traveled  by  the  projectile  base.  The 
other  symbols  were  defined  earlier. 

Finned  projectiles  generally  have  six  fins.  To  reduce  asymmetrical  load.ng 
of  the  fins  as  the  projectile  passes  through  the  vented  region  and  to  limit  cir¬ 
cumferential  pressure  gradients,  brakes  are  usually  designed  with  twelve  columns 
of  holes  so  =  12  will  be  used  in  all  of  the  calculations  below.  The  number 
of  rows  of  holes.  Nr,  will  be  varied  to  obtain  the  vent  area  desired. 
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TABLE  I.  STARTING  SOLUTION  DATA 


120-mm  Cannon 

20-mm  Cannon 

12604.0 

4219.0 

19.59 

0.0878 

15.60 

0.2161 

595.0 

2.545 

5073.6 

3429.0 

1.22 

1.25 

27.2 

30.0 

4.724 

0.7874 

170.9 

56.3 

198.0 

63.7 

0.945 

0.1575 

Parameter 


Pv  (Ibf/in.*) 

C  (Ibm) 

U  ( in. *  ) 

Vpv  (ft/sec) 

y 

T)  (in.Vlb^) 
Ob 

Xg  (in.) 
l-b 


Consider  the  two  hole  patterns  in  Figure  5  where  three  of  the  twelve 
columns  of  each  brake  are  shown;  the  flow  is  from  left  to  right.  The  uniform 
pattern  has  D^.jr,  »  Omax  *  L  *  0.945  inch  for  maximum  efficiency;  to  maintain  the 
same  vent  area,  the  stepped  pattern  used  0^,^  >  0.563  inch,  Dmax  =  1-0  inch,  and 
a  =  0.6  in  the  power  law  expression,  Eq.  (15).  As  indicated  by  the  data  above 
the  patterns,  redistributing  the  area  in  this  manner  increases  the  recoil 
impulse  by  less  than  two  percent.  However,  because  the  total  brake  load  is  sup¬ 
ported  by  a  larger  web  area  in  the  stepped  design,  the  maximum  axial  load  per 
unit  web  area  is  reduced  by  40  percent  (the  maximum  occurs  at  the  instant  the 
projectile  base  is  at  the  muzzle).  The  variation  of  axial  load  per  unit  web 
area  with  distance  along  the  brake  is  shown  in  Figure  6.  The  reduced  stress 
levels  near  the  brake  entrance  should  have  a  favorable  effect  on  fatigue  life. 


UNIFCy?M  PREP  DISTRIBUTION 

RECOIL  IMPULSE  *  3.99S  LB-SEC 

MPXIMUM  PXIPL  LOPO  /  WEB  PREP  =  76.000  PS  I 

oooooooo 

oooooooo 

oooooooo 


STEPPED  PREP  distribution 

recoil  impulse  =  4.050  LB-SEC 

MPXIMUM  PXIPL  LOPO  /  WEB  PREP  =  45. BOO  PS I 

o  o  o  O  O  O  OOOO 
o  o  o  OOOOOOO 
o  o  O  OOOOOOO 

A  comparison  of  a  uniform  hole  pattern  with  a  stepped  pattern. 
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Figure  6.  Axial  load  per  unit  web  area  along  the  brakes 
for  the  two  patterns  shown  in  Figure  5. 


Because  of  the  more  gradual  increase  in  vent  area  with  distance  in  the 
stepped  design,  the  drop  in  propellant  gas  pressure  is  more  moderate,  as  shown 
in  Figure  7.  This  implies  that  the  pressure  gradient  in  the  azimuthal  direction 
will  als  :  be  smaller  so  a  finned  projectile  will  be  exposed  to  a  less  severe 
loading  as  it  passes  through  the  brake. 
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Figure  7.  Pressure  distribution  along  the  brakes  for  the  two  patterns 
shown  in  Figure  5. 

This  example  demonstrates  that  redistributing  a  given  vent  area  along  the 
tube  can  result  in  a  structurally  more  acceptable  design  with  only  a  modest  loss 
in  gasdynamic  performance.  In  practice,  of  course,  the  stepped  pattern  would 
consist  of  a  few  different  hole  sizes  rather  than  continuously  varying  diame¬ 
ters. 
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BRAKE  EFFICIENCY  AND  SCALING 

A  parameter  that  is  useful  for  comparing  the  performance  of  conventional 
muzzle  brakes  is  the  gasdynamic  brake  efficiency,  /3,  defined  as 

0  ~  (^wo  ”  ^w)/(^wo  ~  '"^Pwo)  (27) 

where  is  the  recoil  impulse,  W  is  the  projectile  mass  and  Vp  is  its  muzzle 
velocity.  The  subscript  "wo"  refers  to  the  bare  muzzle  case. 

/3  is  often  found  to  be  remarkably  insensitive  to  the  weapon  caliber  and 
internal  ballistics  and  dependent  only  upon  brake  geometry.  It  is  a  practical 
tool  for  comparing  various  brake  designs  on  a  given  weapon  or  for  estimating 
performance  based  upon  data  obtained  with  a  different  caliber  weapon.  It  will 
be  demonstrated  that  this  utility  carries  over  to  perforated  muzzle  brakes  by 
comparing  the  performance  of  geometrically  similar  brakes  on  a  120-mm  cannon  and 
a  20-mm  cannon.  The  data  in  Table  I  for  the  20-mm  cannon  were  taken  from 
Reference  1. 

The  brakes  have  4,  8,  and  12  rows,  respectively,  of  uniform  diameter  holes 
with  0  *  L.  The  eight-row  pattern  is  shown  in  Figure  5.  The  20-mm  brake  dimen¬ 
sions  were  one-sixth  those  of  the  120-mm  designs. 

In  Figure  8,  brake  efficiency  is  plotted  as  a  function  of  the  vent  area 
ratio,  AR,  which  is  defined  as  the  ratio  of  the  vented  area  to  the  bore  area. 

The  solid  curve  labeled  L/D  =  1  represents  the  calculated  20-mm  data;  the  three 
triangular  symbols  are  the  120-mm  results.  Even  though  the  internal  ballistics 
differs  considerably  for  the  two  cannons,  geometric  similarity  appears  to  be 
sufficient  to  guarantee  nearly  the  same  gasdynamic  performance. 

The  solid  curve  labeled  L/D  =  2  represents  20-mm  calculations  where  the 
wall  thickness,  L,  was  doubled  but  the  number  of  holes  was  kept  the  same.  This 


conforms  to  the  experimental  configurations  of  Reference  1  and  facilitates  com¬ 
parison  with  those  data.  The  asterisk  symbol  is  a  120-mm  calculation  with  L 
fixed,  but  0  halved  to  obtain  L/D  =  2;  the  number  of  holes  was  quadrupled  to 
maintain  the  area  ratio.  This  is  most  likely  how  a  high  aspect  ratio  (L/D) 
would  occur  in  practice.  The  close  agreement  between  the  two  sets  of  calcula¬ 
tions  suggests  that  even  geometrical  similarity  is  not  a  strict  requirement, 
rather,  the  aspect  and  area  ratios  may  be  the  controlling  parameters.  This 
notion  is  reinforced  by  the  remaining  four  symbols  on  the  plot.  These  represent 
a  portion  of  the  20-mm  experimental  data  of  Reference  1  and  correspond  to  brakes 
which  had  only  9  columns  of  holes  and  16,  19,  and  23  rows,  respectively,  rather 
than  12-row  designs  studied  here. 

X  20  MM  experimental  DATA  -  L/0=1.00 
□  20  MM  EXPERIMENTAL  DATA  -  L/D=2.00 
A  120  MM  CALCULATED  DATA  -  L/D=1.00 
„  X  120  MM  CALCULATED  DATA  -  L/Ox2.00 


Figure  8.  Brake  efficiency  as  a  function  of  vent  area  ratio,  AR. 
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Dillon  (ref  2)  recognized  the  relevance  of  these  two  parameters  and  pre¬ 
sented  the  following  formula  for  brake  efficiency: 

3  =  0.208  +  0.312(AR)  -0.145(L/D)  -  0.062(AR)* 

+  0.087(L/D)*  -  0.039(AR) (L/0)  0.006(AR)*  -  0.008(L/D)^  (28) 

The  dashed  lines  in  Figure  8  were  obtained  using  this  expression.  It  overesti¬ 
mates  the  experimental  data  and  the  present  calculations,  particularly  at  small 
area  ratios.  This  happens  because  the  formula  was  based  upon  a  larger  sample  of 
brake  configurations,  some  of  which  were  quite  different  geometrically  from  the 
rest . 


For  example,  one  brake  consisted  simply  of  two  holes  drilled  perpendicular 
to  the  bore  axis  with  diameters  equal  to  the  bore  diameter.  While  this  baffle¬ 
like  configuration  represented  an  interesting  comparison  for  the  more 
"conventional"  perforated  brakes,  its  inclusion  in  the  formula  accounts  for  the 
two  dashed  lines  converging  to  a  high  /$  at  AR  =  2  in  the  plot. 

Another  configuration  had  vents  only  on  the  sides  of  the  tube.  The  lateral 
venting  arrangement  was  an  attempt  to  reduce  possible  obscuration  effects  asso¬ 
ciated  with  axisymmetric  venting.  As  noted  in  the  Introduction,  full-scale 
testing  of  the  idea  resulted  in  damage  to  finned  projectiles.  It  did  produce 
somewhat  higher  efficiencies  in  the  20-mm  tests,  however,  and  this  contributes 
to  the  high  estimates  of  the  formula. 

Nevertheless,  Dillon's  idea  of  possibly  representing  a  wide  range  of  per¬ 
forated  configurations  with  a  single  formula  has  obvious  appeal  and  can  be  more 


This  was  obtained  using  the  20-mm  calculations  of  Figure  8.  All  of  the  data  in 
Reference  1  for  axisymmetric  configurations  with  aspect  ratios  between  1  and  2 
are  shown  as  well  as  the  uniform  and  stepped  120-mm  data  calculated  here.  This 
information  is  summarized  in  Tables  11  and  III. 

❖  20  MM  experimental  DATA  -  L/D=0.9S 

X  20  MM  experimental  DATA  -  L/0=I.0O 

O  20  MM  EXPERIMENTAL  [DATA  -  L/D=l.26 

Z  20  MM  EXPERIMENTAL  DATA  -  L/Oal.89 

□  20  MM  EXPERIMENTAL  DATA  -  L/D=2.00 

A  120  MM  CALCULATED  DATA  -  L/D=1.00 

+  120  MM  CALCULATED  DATA  -  L/0=1.09 

_  »  120  MM  CALCULATED  DATA  -  L/D=2.00 


Figure  9.  Correlation  of  computed  and  experimental  brake  efficiencies. 

The  model  numbers  in  Table  II  correspond  to  those  in  Reference  1.  However, 
the  values  of  0  are  slightly  lower  than  reported  there.  In  computing  the  brake 
efficiency  from  the  impulse  measurements,  Dillon  used  the  muzzle  velocity 
obtained  with  the  brake  in  place  rather  than  the  bare  muzzle  velocity.  To  be 
consistent  when  comparing  the  various  brake  configurations,  the  latter  was  used 
here. 


20 


TABLE  II.  SUMMARY  OF  EXPERIMENTAL  DATA  FROM  REFERENCE  1 


Model 

Pattern 

AR 

L/D 

■m 

1 

staggered 

6.32 

1.26 

70.4 

2 

staggered 

5.00 

1.89 

51.3 

3 

uniform 

4.13 

1.00 

56.8 

4 

uniform 

4.13 

2.00 

44.3 

6 

uniform 

4.89 

2.00 

51.1 

7 

uni  form 

5.99 

2.00 

54.4 

8 

staggered 
&  stepped 

3.38 

0.95 

50.4 

TABLE  III.  SUMMARY  OF  120-mm  CALCULATIONS 


Model 

Pattern 

AR 

L/D 

1 

1 

A 

uni  form 

1.92 

1.00 

32.8 

B 

uniform 

3.84 

1.00 

52.6 

C 

uniform 

5.76 

1.00 

67.9 

0 

stepped 

3.82 

1.09 

50.5 

E 

stepped 

5.74 

1.09 

65.0 

F 

uniform 

1.92 

2.00 

24.6 

The  staggered  patterns  in  Table  II  had  every  other  row  of  12  holes  rotated 
15  degrees  with  respect  to  the  adjacent  rows.  The  uniform  pattern  refers  to  the 
column  arrangement  shown  in  Figure  5.  Staggering  is  expected  to  reduce  azi¬ 
muthal  pressure  gradients  near  projectile  fins.  It  does  not  appear  to  affect 
gasdynamic  performance,  but  calculations  show  that  the  column  pattern  produces 
somewhat  lower  stresses  for  a  given  load  (ref  lA). 


The  stepped  pattern  in  Table  II  had  three  rows  of  holes  with  l/16-inch  diame¬ 
ters,  three  rows  with  1/8-inch  diameters,  and  four  rows  with  3/16-inch  diame¬ 
ters.  The  center-to-center  spacing  in  the  axial  direction  was  constant,  as  in 
the  stepped  pattern  of  Figure  5,  although  the  rows  were  staggered.  Each  row  had 
12  holes. 

To  calculate  an  average  hole  diameter  for  the  stepped  patterns,  the 
following  formula  was  used: 

N  N 

0  =  I  MnAn  /  }  NnAp  (30) 

n=l  n=l 

where  N  is  the  number  of  different  hole  sizes  and  Np,  is  the  number  of  rows  with 
diameter  0^  and  area  Ap.  Each  diameter  is  weighted  by  its  area  in  the  belief 
that  a  large  hole  contributes  more  to  brake  performance  than  a  small  hole.  The 
weighting  factor  shifts  the  data  closer  to  the  curve  in  Figure  9. 

The  close  correlation  of  the  tabulated  data  by  the  predicted  curve  supports 
Dillon's  observation  that  brake  performance  is  primarily  dependent  on  vent  area 
ratio  and  aspect  ratio.  Stepped  patterns  are  included  by  using  a  suitably 
weighted  average  diameter  to  compute  the  aspect  ratio. 

CONCLUSIONS 

The  influence  on  brake  performance  of  including  covolume  in  the  three- 
dimensional  calculations  was  found  to  be  too  small  to  warrant  the  added  com¬ 
putation.  The  covolume  term  should  be  included  in  the  one-dimensional  code, 
however,  to  be  consistent  with  the  internal  ballistics  code  used  to  generate  the 
starting  solution. 


Redistributing  a  given  vent  area  along  the  tube  using  a  set  of  holes  with 
varying  diameters  can  result  in  a  structurally  more  acceptable  design  with  only 
a  modest  loss  in  gasdynamic  performance.  It  may  also  present  a  more  favorable 
environment  for  finned  projectiles  by  reducing  axial  and  azimuthal  pressure  gra¬ 
dients  in  the  vented  portion  of  the  tube. 

Brake  performance  is  primarily  dependent  on  vent  area  ratio  and  aspect 
ratio.  The  dependence  can  be  expressed  by  a  simple  formula  for  axisymmetric 
designs.  Stepped  patterns  are  included  by  using  a  suitably  weighted  average 
diameter  to  compute  the  aspect  ratio. 
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